Is the Quantum Zeno Effect Evolution's Choice for the Avian Compass? 
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Magnetic-sensitive radical-ion-pair reactions are understood to underlie the biochemical magnetic 
compass used by avian species for navigation. Radical-ion-pair reactions were recently shown to 
manifest a host of quantum-information-science effects, like quantum jumps and the quantum Zeno 
effect. We here show that the quantum Zeno effect immunizes the magnetic and angular sensitivity 
of the avian compass mechanism against the deleterious and molecule-specific exchange and dipolar 
interactions. 

PACS numbers: 82.20.-w 



"In the history of natural selection, did nature ever 
come across a way to use quantum weirdness?" This is 
a question claimed [l| to have an affirmative answer, at 
least in regard with the apparent ability of photosyn- 
thetic antennae to efficiently guide the excitonic energy to 
the photosynthetic reaction center. The possibility of bi- 
ological processes exhibiting non-trivial quantum effects, 
ordinarily thought to be suppressed in the decoherence- 
prone biological environment 0], has attracted an in- 
creasing attention in recent years. For example, signif- 
icant experimental Q and theoretical [1] progress has 
been recently made on elucidating the role of quantum 
coherence in the workings of photosynthetic antennae 
cornplexes. In a different front, radical-ion-pair reactions 
[^, 0] have been recently shown [7] to exhibit the full 
machinery of concepts and physical effects familiar from 
quantum information science. Radical-ion pairs play a 
fundamental role in a series of biologically relevant chem- 
ical reactions [1, [§] , and are understood to underlie the 
biochemical magnetic compass use d by avian species to 
navigate in earth's magnetic field [lO|, as corroborated 
by several recent experiments jllj. 

Radical-ion pairs are formed by a charge transfer fol- 
lowing a photoexcitation of a donor-acceptor dyad, lead- 
ing to two molecular ions and two unpaired electrons. 
The latter can either be in the spin singlet or in the 
spin triplet state. Magnetic interactions with the ex- 
ternal magnetic field and hyperfine interactions with 
the molecule's magnetic nuclei bring about a coherent 
singlet-triplet oscillation. At some random instant in 
time the reaction is terminated, since the radical-ion-pair 
undergoes charge recombination, leading to the reaction 
products. Spin conservation enforces spin selectivity of 
the recombination process, i.e. singlet (triplet) radical- 
ion pairs recombine to singlet (triplet) neutral products. 
Moreover, anisotropic hyperfine interactions within the 
molecule render the reaction yields dependent on the in- 
clination of the external magnetic field with respect to a 
molecule-fixed coordinate frame. The avian magnetore- 
ception mechanism exactly rests on both aforementioned 
aspects, i.e. the reaction yield's magnetic and angular 



sensitivity. 

Interestingly, intra-molecule magnetic interactions are 
more complicated. Both spin-exchange and long-range 
dipolar interactions affect reaction dynamics, to an ex- 
tent dependent on the particular molecular structure. In 
this respect, it has been recently shown [l^] that the pres- 
ence of exchange and/or dipolar interactions significantly 
suppresses the magnetic and angular sensitivity of the 
reaction yields, thus severely degrading the mechanism's 
functionality. Along the same lines, it was concluded 
[l3 | that only when the molecular parameters determin- 
ing J and D, the exchange and dipolar couplings, are 
fine-tuned so that the effects of these two interactions 
cancel each other, is the magnetic and angular sensitivity 
of the reaction restored. Although a remarkable result, 
it is hard to imagine that Nature has conjured up such 
a fortuitous cancellation for a functionally important bi- 
ological sensor. Even more so since J and D depend 
exponentially on molecule-specific parameters, like the 
donor-acceptor distance. 

However, it was recently shown It'I that the fundamen- 
tal description of radical-ion-pair reactions following from 
quantum measurement theory reveals the actual presence 
of non-trivial quantum effects, such as the quantum Zeno 
effect [l^. We will here show that when the quantum 
Zeno effect is manifested, the reaction's magnetic and 
angular sensitivity is practically independent of the pres- 
ence or not of exchange and/or dipolar interactions. This 
realization has profound implications for the robustness 
of this biological sensor, i.e. a non-trivial quantum ef- 
fect renders the sensor insensitive to molecule-specific pa- 
rameters [l2|, such as donor-acceptor distance (affecting 
the exchange coupling and the long-range dipolar cou- 
pling) and the inter-radical medium and the particular 
electronic structure (affecting the exchange coupling). 

The spin state of the radical-ion-pair is described by a 
4rt-dimensional density matrix p , where the factors 4 and 
n are the spin multiplicities of the two electrons and the 
molecule's nuclear spins, respectively. The unconditional 
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evolution of p is given f?! by the Liouville equation 

dp/dt = ~i[Hm, p] - {ks + kT){Qsp + pQs - '^QspQs) 

(1) 

where Qs and Qs are the singlet and triplet state pro- 
jection operators, respectively. The first term in (1) is 
the unitary evolution due to the magnetic interactions, 
whereas the second term describes the measurement- 
induced evolution due to the singlet (rate ks) and triplet 
(rate fcs) recombination channel. The reaction is ter- 
minated at a random instant in time according to the 
quantum jump equations [lij 



dps = 2ks{Qs)dt 
dpT = 2kT{QT)dt 



(2) 



which give the probability of singlet and triplet recom- 
bination taking place within the time interval between t 
and t + dt. 

The magnetic Hamiltonian Hm — Hz + Ti-hf + 'Hex 
is composed of Hz, the Zeeman interaction of the two 
unpaired electrons (nuclear Zeeman interaction is neg- 
ligible) with the external magnetic field, Huf, the hy- 
perfine couplings of the electrons with the surrounding 
nuclear spins, and finally the spin-exchange interaction, 
Hex- For the transparency of the following discussion 
we will ignore the dipolar interaction as its inclusion 
leads to exactly the same conclusions. The Zeeman in- 
teraction Hamiltonian that will be used for the study 
of the magnetic sensitivity is Hz,magn — ^{siz + S22), 
where the magnetic field of magnitude B is assumed to 
be in the x-axis (w = ^B, with 7=1.4 MHz/G). For 
the study of the angular sensitivity we take the mag- 
netic field, again of magnitude B, to be in the x-y plane, 
hence the Zeeman interaction term will be Hz,ang = 
ujcos(j){six + S2x) + ujsm(j>{siy + S2y)- In the follow- 
ing we will consider the simplest physically realizable 
radical-ion-pair containing just one spin-1/2 nucleus (in 
which case dim(p)=8), hence the hyperfine interaction 
Hamiltonian is Hhf = I • A • Si, where A is the hyper- 
fine coupling tensor of the single nuclear spin I existing 
in e.g. the donor molecule with the donor's unpaired 
electron. We will consider the simplest case where the 
hyperfine tensor is diagonal with one non-zero compo- 
nent Axx = a, to provide for the angular sensitivity on 
the X — y plane. Thus Hhf = asixlx- Finally the spin- 
exchange Hamiltonian is Hex = Jsi -82. In the following 
we will calculate the magnetic and angular sensitivity of 
the reaction for two regimes: (i) the "traditional" regime 
where both recombination rates fcg and fc-r are equal and 
small compared to the magnetic frequency scale h (set 
by a and lu), i.e. ks — kx < a,LU. It is in this regime 
that almost all calculations have been performed based 
on the previous, phenomenological theory We then 
study the regime (ii) where kr S> ks and kr ^ uj, i.e. 
when the quantum Zeno effect is manifested, since the 
total measurement rate k = ks + kx is dominated by 



the triplet recombination channel, strongly re-projecting 
the radical-ion-pair's spin state to the singlet state, from 
which it is considered to start out, i.e. the initial condi- 
tion is p{t = 0) = Qs/2 (which satisfies Tr{p} = 1 and 
Tr{pQs} = 1). We will show that in the regime (i) the in- 
clusion of the spin exchange interaction indeed degrades 
the magnetic and angular sensitivity of the reaction, as 
has already been analyzed However, regime (ii) ex- 
hibits a significant magnetic and angular sensitivity with 
their dependence on the exchange coupling J being al- 
most entirely suppressed. 

In Fig. 1 (insets) we show the triplet yield Yt of the 
reaction as a function of the external magnetic field (us- 
ing Hz,magn) and the field's angle (using Hzmtuj)- The 
triplet yield is calculated from 

100 f°° 

^T[%] = ^y^ dPTmit) (3) 

where Nq is the initial number of radical-ion pairs (start- 
ing out in the singlet state) and N{t) the remaining 
number of radical-ion pairs at time t, evolving as [7[ 
dN = —N{dps + dpx)- That is, in the time interval be- 
tween t and t + dt the number of created triplet products 
will be the number of radical-ion-pairs existing up to that 
time, N{t), times the probability dpx for triplet-channel 
recombination. This is summed in ([3]) until all molecules 
have recombined one way or another. The probabilities 
dps and dpx are calculated from 121), whereas the expec- 
tation values {Qs) and (Qt) entering ^ are obtained by 
evolving the density matrix according to ([T]) . The mag- 
netic sensitivity of the reaction at earth's field of interest 
for the avian compass is just the slope dYr/dB calcu- 
lated at 5=0.5 G. The smallest measurable change of 
the magnetic field (magnetic precision), 6B, and the pre- 
cision Scf) of the determination of the magnetic field angle 
with respect to the molecule's x-axis (angular precision 
or heading error) both follow from the previous calcu- 
lations if the smallest measurable reaction yield change, 
SYt, is known. It thus follows that 



SB = 



SYt 



\dYT/dB\B=0.5 G 

6Yt 



(4) 
(5) 



(Yt .max Yt .mm 

where Ir.mm and YT^max are the minimum and maxi- 
mum values of the yield Yt {(f) and = 90° is the an- 
gular width of the full swing between YT^min and YT.max ■ 
To make further progress towards the absolute mag- 
netic and angular sensitivity of the avian compass, the 
smallest "measurable" change in the reaction yield 6Yt 
must be known or estimated. Obviously SYt depends 
on the particular realization of the biochemical mecha- 
nism transducing the radical-ion-pair reaction yield to a 
physiological signal. On rather general grounds it has 
been shown [15| that SYt is connected to Nu, the num- 
ber of neuronal receptors sensitive to the radical-ion-pair 
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FIG. 1: Magnetic and angular precision. For all calculations 
the hyperfine coupling is a =5 G. (a) Inset shows an exam- 
ple of the magnetic sensitivity of the triplet yield plotted in 
the quantum Zeno regime for two values of the exchange cou- 
pling J. For higher (lower) values of the exchange coupling 
the plots are similar but appear at a lower (higher) offset. 
The yellow band is the 0.05% error band corresponding to 
the measurement precision of the reaction yield. Main graph 
shows the magnetic precision, SB, extracted from the slope at 
B =0.5 G of plots like the ones in the inset, as a function of 
the exchange coupling J for the two regimes of the recombina- 
tion rates, (b) Similar calculation as in (a), but using Tiz.ang 
instead of Ti.z,magn, in order to derive the angular precision 
or heading error. 



reaction product molecules, by (SYt)^ — We 
choose iVfl — 1.6 x 10'', in order to set SYt at the value 
SYt — 0.05%. The chosen value of Nr and hence SYt 
is realistic [l5| and has the consequence that it sets the 
magnetic precision at zero exchange coupling (shown in 
Fig. la) at the value ofSB « 0.01 G, i.e. at 2% of earth's 
field. This level of magnetic precision is understood [l3| 
to be actually the case for several avian species. It is 
stressed, however, that the following considerations are 
qualitatively independent of the particular value of SYt, 
which just sets the absolute scale of the derived magnetic 
and angular precision. It is clearly seen in Fig. la that in 
the traditional regime (i) the magnetic precision plunges 
from SB « 0.01 G at J = to over 1 G at J=15 n G. 
Even at J— 8 G the change of earth's field is barely mea- 
surable. On the contrary, in the quantum Zeno regime 
(ii) the magnetic precision drops only to SB « 0.07 G at 



J=15 G. Similarly, as shown in Fig. lb, the angular pre- 
cision in regime (i) drops dramatically with increasing J, 
with a complete loss of heading information for J=10 G 
or higher. In contrast, in the quantum Zeno regime the 
angular precision of about S(j) — 20° at the highest value 
of the exchange couplin g is actually at the level of ex- 
perimental observations [iq of the heading error of the 
avian compass. At the same time the heading error is 
practically unaffected by the presence or not of exchange 
and dipolar interactions, varying between Scj) = 3.5° and 
S(f) — 20.5° in the full span of the exchange coupling J. 

We will now explain the robust magnetic and angu- 
lar precision in the quantum Zeno regime. This follows 
by considering the behavior of the eigenvalues of the 
master equation ([T]), which are obtained by diagonaliz- 
ing the matrix Ai (of dimension dim(p)^) that satisfies 
dp/dt = Mp, where p is a column matrix containing all 
matrix elements of p . The resulting eigenvalues are of 
the form —A + iil, with A > being the decay rate and 
f2 the oscillation frequency of the particular eigenmode. 
As is in general the case with the quantum Zeno effect 



17| . some of the eigenvalues have decay rates increasing 
with the measurement rate fc as A ~ /c, while the others 
(responsible for the quantum Zeno effect) decrease with 
k as Xqz ~ /i^/fc, where h is the characteristic frequency 
scale of the system, here determined by the magnetic 
Hamiltonian Tim- 

Now, the exchange Hamiltonian can be written (up to 
an additive constant) as Hex — —JQs- Furthermore, as 
is known from quantum measurement theory [l8| , the de- 
terministic evolution of the system's quantum state (due 
to the unitary Hamiltonian evolution and the measure- 
ment of Qs with rate k ^ ks + kT) is generated by the 
non-hermitian operator IC = Tim — ikQs- It is easily seen 
that if Ti-j^^ is the magnetic Hamiltonian without the ex- 
change interaction Tiex^ then /C = H.^^^ — i{k — iJ)Qs, 
i.e. the inclusion of the exchange interaction is equiv- 
alent to replacing k with with an imaginary measure- 
ment rate k — iJ. We can now complete the argument 
as follows: the eigenvalues with a real part that scales 
as A ~ fc pick up an oscillation frequency (in addi- 
tion to n) of — J, the effect of which roughly averages 
out. On the other hand, the eigenvalues with the quan- 
tum Zeno scaling Xqz ~ h'^/k suffer a change in their 
real part which becomes (since in our case k/J <^ 1) 
^'qz ~ h'^k/J^ <^ Xqz- Thus the spin state evolution is 
slowed down, the radical-ion-pair is projected stronger 
to the singlet state leading to an increased singlet yield, 
or equivalently a decreased triplet yield, which is ex- 
actly what is observed in Fig. 1 (insets). The above 
argument is pictured in Fig. 2a, which shows the time 
evolution of N(t), the number of existing (non-reacted) 
radical-ion pairs. The reaction is considered to be ter- 



minated when N{t)/N{0) 



5 X 10- 



i.e. when the 



reaction yield is known to within SYt- It is clearly seen 
that in the quantum Zeno regime, the reaction time de- 
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FIG. 2: Reaction Dynamics, (a) Time evolution of the num- 
ber of radical ion pairs for the two regimes of the recombi- 
nation rates, for different values of J. The dashed line at 
N = 5 X 10^ signifies the "termination" of the reaction, i.e. 
the point when the remaining radical-ion pairs are 0.05% of 
the initial number N{t = 0) = 10^°. (b) Time evolution of 
{Qs(t)> and (QtW) for 5=0.5 G and J=10 G. (c) Triplet 
yield as a function of the magnetic field for J=10 G. 

pends on J in the way outlined before. On the contrary, 
when ks = kx = k, the change in the number of radical- 
ion-pairs during the time interval dt is (Qs + Qt = 1) 
dN = -Ndt{2ks{Qs) + 2kT{QT)) = -2KdtN, hence 
N{t) = iV(0)e~^''*, i.e. the reaction time is proportional 
to 1/k and independent of J. The result is that dur- 
ing the short reaction time in the traditional regime, the 
triplet probability Qt has not increased appreciably (Fig. 
2b), and the triplet yield is small, as shown in Fig. 2c. 
In contrast, in the quantum Zeno regime the reaction has 
enough time to "sample" large values of Qt and lead to 
a triplet yield about an order of magnitude higher, hence 
the higher sensitivity in this regime. In other words, as 
seen in Fig. 2c, the relative change SYt/Yt of the triplet 
yield with the magnetic field is roughly the same in both 



cases, but the absolute value of It differs by a factor 
of 10, leading to respectively high slopes dYT/dB and 
dYT/d(j). To summarize, the quantum measurement dy- 
namics inherent in the recombination process of radical- 
ion pairs result in " delocalization" of the electron spin 
state at long times, as evidenced in Fig. 2b. The asym- 
metric {kT S> ks) recombination rates result in the ac- 
tual manifestation of the quantum Zeno effect and the 
J-dependence of the reaction dynamics. The interplay 
of these two effects provides for the robust magnetic and 
angular sensitivity in the presence of the exchange inter- 
action. 

Whether Nature has eng ineered molecules realizing a 
fine-tuned cancellation [l^] of the adverse effects of ex- 
change and dipolar interactions, or on the other hand, has 
evolved the avian compass into operating at the quantum 
Zeno regime, remains to be discovered. 
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